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Abstract 

An effective monopole action for various extended monopoles is derived from 
vacuum configurations after abelian projection in the maximally abelian gauge 
in SU (2) QCD. The action appears to be independent of the lattice volume. 
Moreover it seems to depend only on the physical lattice spacing of the renor- 
malized lattice, not on (5. Entropy dominance over energy of monopole loops 
is seen on the renormalized lattice with the spacing 6 > 6c — 5.2 x 10^^ A^^. 
This suggests that monopole condensation always (for all jS) occurs in the 
infinite-volume limit of lattice QCD. 
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To clarify color confinement mechanism is still important in particle physics and the 
effect dual to the Meissner effect is believed to be the mechanism This picture is 

realized in the confinement phase of lattice compact QED [|]-|^. 

To apply the idea to QCD, we have to find a color magnetic quantity inside QCD. The 
'tHooft idea of abelian projection of QCD is very interesting. The abelian projection of 
QCD is to extract an abelian gauge theory by performing a partial gauge- fixing such that 
the maximal abelian subgroup remains unbroken. Then QCD can be regarded as an abelian 
gauge theory with magnetic monopoles and electric charges, 't Hooft conjectured that the 
condensation of the abelian monopoles is the confinement mechanism in QCD 0. 

There are, however, infinite ways of extracting such an abelian theory out of QCD. It 
seems important to find a good gauge in which the conjecture is seen to be realized in a 
coupling region where reliable calculations can be done at present. A gauge called maximally 
abelian (MA) gauge has been shown to be very interesting 0-3. In the MA gauge, there are 



phenomena which may be called abelian dominance PJTO||. Moreover the monopole currents 
kf^{s) defined similarly as in compact QED ^ are seen in the MA gauge to be dense and 
dynamical in the confinement phase, while dilute and static in the deconfinement phase . 
Especially the string tension can be reproduced by the monopole contributions alone ||TT],|T2 . 

Although these data support the 'tHooft conjecture, there is not a direct evidence of 
abelian monopole condensation in the QCD vacuum. In the case of compact QED, the 
exact dual transformation can be done and leads us to an action describing a monopole 
Coulomb gas, when one adopt the partition function of the Villain form PJT^-[15| . Monopole 



condensation is shown to occur in the confinement phase from energy-entropy balance of 
monopole loops. 

In the case of QCD, however, we encounter a difficulty in performing the exact dual 
transformation. Hence we tried to do it numerically first by fixing an effective U(l) action. 
Namely, the abelian dominance suggests that a set of U{1) invariant operators are enough 
to describe confinement after the abelian projection. Then there must exist an effective 
U{1) action describing confinement. As reported in [^], however, to derive a compact form 



of such a U{1) action in terms of abelian link variables was unsuccessful in the scaling 
region. Here we propose a method to directly determine a monopole action numerically 
from the vacuum configurations of monopole currents given by Monte-Carlo simulations in 
the MA gauge |T6| , p^ . This can be done by extending the Swendsen method ||18| which was 
developed in the studies of the Monte-Carlo renormalization group. We have shown recently 
that the method is successfully applied to compact QED |T^. When applied to SU(2) QCD 
in the MA gauge, it is suggested that entropy dominance over energy and then monopole 
condensation always occur in the infinite- volume limit of lattice QCD. A preliminary report 
was published in ||l6| , p!7| . 

The MA gauge is given on a lattice by performing a local gauge transformation 



U{s, fi) = V{s)U{s, fi)V~\s + fi) 



(1) 



such that a quantity 



R = Y1 Trfdst/ (s, iJ,)(TsU\s, /i) 



(2) 



is maximized [0. In the continuum limit, this corresponds to a gauge {d^+igA^^){Aj^—iA'j^) = 
0. Such a gauge fixing keeps a residual U(l) invariance. 

After the gauge fixing is done, an abelian link gauge field is extracted from the SU{2) 
link variables as follows; 

( (l-|c(.,/i)P)V2 



-c*(s,/x) 



V 



V 







(3) 



C(s,/i) (l-|c(s,/i)|2)V2 

where u{s,fi) = exp{i9{s, fi)) represents the abelian link field and c{s,fi) corresponds to 
charged matter fields. 

A plaquette variable 6^y{s) E (— 47r,47r) is given by a link angle 9{s,fi) E [— vr, tt) as 
(^fj,u{s) = 9{s,n) + 0{s + fi^u) — 9{s + z>, /i) — 9{s, z/). The plaquette variable 0^^{s) can be 
decomposed into 



(4) 



where 6^y{s) E [— vr, vr) is interpreted as the electromagnetic flux through the plaquette and 
^nu{s) can be regarded as a number of the Dirac string penetrating the plaquette. DeGrand 
and Toussaint defined a monopole current as 

^/^(■s) = ]j(^^lupcT^ump„{s + /i), (5) 

where d is the forward difference. It is defined on a link of the dual lattice (the lattice with 
origin shifted by half a lattice distance all in four directions). k^j_{s) can take only the values 
0, ±1, ±2 and satisfies 

where d' is the backward derivative on the dual lattice. 

Long-distance behaviors are expected to be important in the confinement phase of QCD. 
Hence we consider extended monopoles of the type-2 as shown in Fig. 0. |^ The 
extended monopole of the type-2 has a total magnetic charge inside the cube and is 
defined on a sublattice with the spacing b = na, a being the spacing of the original lattice. 

kP{s) = le,^,.d,m(;J{s + fi) (7) 

n-l 

= X! k^{ns + {n - l)fl + iu + jp + la), (8) 

i,j,l=0 
n-l 

"^ia^s) = J2 ^pAns + ip + ja). (9) 

i,j=0 

The definition of the type-2 extended monopoles corresponds to making a block spin trans- 
formation of the monopole currents with the scale factor n. We call the sublattice as a renor- 
malized lattice. We derived the coupling constants for 1^, 2^, 3^, and 4'^ extended monopoles. 
A theory of monopole loops is given in general by the following partition function 

oo 

Z = il[ E )iI[So>,k,is),o)eM-S[k]), (10) 

fc,i(s)=-oo s 



^ Type-1 extended monopoles are not adopted here, since the separation of the Dirac string 
is problematic. 
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where k^{s) is the conserved integer- valued monopole current defined above in (||) or and 
S[k] is a monopole action describing the theory. Consider a set of all independent operators 
which are summed up over the whole lattice. We denote each operator as Si[k]. Then the 
action can be written as a linear combination of these operators: 

s[k] = Y.mk], (11) 

i 

where fi are coupling constants. 

Let us now determine the monopole action using the monopole current ensemble which 
are calculated from vacuum configurations generated by Monte-Carlo simulations. Since the 
dynamical variables here are kf^{s) satisfying the conservation rule, it is necessary to extend 
the original Swendsen method Consider a plaquette (s',/i', u') instead of a link on the 
dual lattice. Introducing a new set of coupling constants {fi}, define 

^ r^j ^ J2m=-oo Si[k'] exp(— J2i fiSi[k']) 

EM=-ooeM-J:^fiSi[k']) ' 

where k'^{s) = k^{s) + M((5s,s'5^^' + 5s,s'+[l'5^,u' - ^s.^'+f-'^^M' ~ ^s,s'5^v')- When all fi are equal 
to fi, one can prove an equality {S.-^ = (Si), where the expectation values are taken over the 
above original action with the coupling constants {/i}: 



(n\k]\ - (n.,MEg(.)=-oc)(ns^a;fc^(s),o)Q[fc]exp(-E»/»'^n--./ . . 

When there are some fi not equal to fi, one may expand the difference as follows 

{s, - s,) = J2(s;s] - s.s,){f, - /,), (14) 

where only the first order terms are written down. This allows an iteration scheme for 
determination of the unknown constants /j. For more details, see ||T^ where ( [T^ ) is derived 
explicitly and is applied to compact QED. 

Practically we have to restrict the number of interaction terms. We adopted 12 types 



^ All possible types of interactions are not independent, since d'^k^{s) = 0. We can get rid of 
almost all interactions between different components of the currents from the quadratic action by 
use of the conservation rule. 
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of quadratic interactions listed in Table | in most of these studies. To check the validity of 
the truncation, we also included more quadratic and even quartic interactions in some cases, 
but found little difference. 

Monte-Carlo simulations were done on 8^, 12^, 14^, 16^, 18^, 20^, and 24^ lattices for < 
(3 < 3.0. 

Our results are summarized in the following. 

1. The coupling constants /i ~ /s for 1^ monopoles are plotted (white symbol) for small 
(3 in Fig. |^. As shown in |^, an effective U{1) action is determined for small (3. It 
takes a Wilson f/(l) action with the effective coupling constant [3^ = (3/2. In Fig. |^, 
we have also plotted the data (black symbol) obtained in the case of the Wilson f/(l) 



action |19] with the coupling constant substitution. Equivalence of the two theories 



for small j3 is seen also in the form of the monopole action. 

2. For (3 > 2.4, extended monopoles are considered. ^ The coupling constants fi are fixed 
for not so large value of b ,i.e., b < 1.3 x 10^^ (A^,)^^. The first important observation is 
that the coupling constants determined are almost independent of the lattice volume 
as seen from Fig. ^ We see /i is dominant and the coupling constants decrease rapidly 
as the distance between the two monopole currents increases. 

3. Since the monopole action is determined, let us consider energy-entropy balance of 
monopole loops as done in compact QED [Q. Monopole world lines form closed loops 
on the dual lattice. Define nilk] as the number of monopole loops of length L in a 
monopole configuration {k^{s)}. The contribution of monopole loops of length L is 
represented by the average of niik] : 



E{fc} nL[k]e 



^S[k] 



(Mk]) = > (15) 



^ In the case of compact QED, a monopole action could be determined only for 1^ monopoles. 
For extended monopoles, the iteration to fix the action was not convergent. 
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where 



E=(n E mh^.is),o)- (16) 

Monopole loops of length L take various forms. Choose a monopole loop Cl of length 
L and define 

[ 1 iiCie {k} 
ScAk] = { (17) 
otherwise . 



nL[k] is written as 



nL[k]=j:6cA% (18) 

Cl 



where J2cl denotes the sum over all loops of length L. Then the average (15) is 



{nL[k]) = j:{Scdk]). (19) 

Cl 

(5(7^ [fc]) depends not only on the loop length L but also on the shape of Cl- 

One may define the 'energy' E{L) of monopole loops of length L as the average of 
(5(7^ [/c]) over all loops of length L: 

^-E{L) ^ HcA^cA^) ^20) 
1 

/X r,n_ i:w^c.[fc]e-^W 

(^c.[A;])-^-— ^j^. (21) 

Then we get 

K[fc])=p(L)e-^(^), (22) 

where p(L) = X^c^ 1 is the total number of monopole loops of length L on a dual 
lattice and Inp(L) may be regarded as the 'entropy' of monopole loops of length L. 

Now we write 
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where S{Cl) is the contribution of a loop Cl alone to the monopole action. For large L, 
monopole currents in a loop are distributed randomly in average. If the action consists 
of quadratic terms alone as assumed in our analysis, interaction terms between different 
currents would cancel. S{Cl) may be approximated by the self-energy part: 

S{Cl) ~ /i X L, (24) 

where we assumed the dominance of currents with a unit charge = 1. 

On the otherhand, Y.{k} Sci^[k]e~^^''^'^^^'-'^^ has two types of contributions, i.e., l)those 
from other loops in the presence of Cl and 2)the interactions between Cl and other 
loops. The latter may be cancelled again from randomness when L is large. When we 
fix Cl and make the lattice volume larger, the ratio of the number of links occupied 
hj Cl to total link number becomes smaller. Hence the Cl dependence of the former 
may be negligible in the infinite- volume limit. We may approximate for large L 

(5c,W)oce-^^^ (25) 

which leads to 

E{L) ~ f\L + const. (26) 

If we consider non-backtracking random walks on a four dimensional dual lattice, the 
number of loops of length L through a given site behaves asymptotically as 7^ for large 
L. We get 



We have made a histogram analysis of the magnetic charge distribution of each extended 
monopole. The case with |fc^(s)| = 2 is much (less than 5%) suppressed in comparison with 
the |A;^(s)| = 1 case and charges with |A;^(s)| > 3 barely appear even for 4^ extended monopoles in 
the confinement phase. 
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p(L) oc 7^. 



(27) 



From (1^) and (pT]) , we obtain 

{nL[k]) oc e^'-^'-^i)^ (28) 

for large L. If /i < ln7, the 'entropy' dominates over the 'energy' and the long 
monopole loops become dominant. Namely the monopole condensation occurs. The 
same thing happens also in compact QED 

We plot fi versus (5 for various extended monopoles on 24'^ lattice in comparison with 
the entropy value ln7 for the infinite volume in Fig. ^. Each extended monopole has 
its own 13 region where the condition /i < ln7 is satisfied. When the extendedness is 
bigger, larger (3 is included in such a region. 

4. The behaviors of the coupling constants are different for different extended monopoles. 
But if we plot them versus h , we get a unique curve as in Figs. | and |^. The coupling 
constants seem to depend only on 6, not on the extendedness nor (3. This suggests 
the existence of the continuum limit and the monopole action in the continuum may 
be similar to that given here. From Fig. we see the self-energy per monopole loop 
length fi on the renormalized lattice decreases as the length scale h increases. A critical 
length hf. ~ 5.2 x exists at which the /i value crosses the ln7 line. 

Together with the first observation of the volume independence, we may get the fol- 
lowing important conclusion in the infinite- volume limit. One can take always (for 
all (3) a renormalized lattice with a spacing 6 > 6c by making a block spin trans- 
formation of the monopole currents with a sufficiently large scale factor. Hence 
the QCD vacuum in the infinite-volume limit of lattice QCD is always (for all (3) 
in the monopole condensed phase. This is the most important observation of this re- 
port. 
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5. What is the meaning of the critical length 6c? When the lattice distance b of the 
renormalized lattice is larger than the critical length b^, the condensation of monopole 
loops is shown from the above 'energy-entropy' balance. If monopoles are physical 
quantities, it must have a finite energy and a finite physical size Rm- Q When the 
lattice distance b is less than the physical monopole size Rm, one can not regard 
monopoles as a point particle running a world line on such a renormalized lattice. 
Hence the 'energy' and the 'entropy' of such monopoles can not be represented as in 
( p6l) and (|27|) . The condition fi < ln7 can not be applicable to such monopoles. The 
above consideration suggests R 0{b,). 

6. Using the data of fi, we try to fix the explicit form of the monopole action. Although 
the coupling constants for largely separated currents are too small with large errors, 
the monopole action may be fitted by the action predicted theoretically by Smit and 

sijs ig: 

1 Att 

S[k] = Y.m,bk,{s)k,{s) + -{-^fY.Us)D{s - s')k,{s'), (29) 
where g{b) is the SU(2) running coupling constant 



The scale parameter determined here is A ~ 42A/^ which is different from that in [|15 



D{s) is a modified lattice Coulomb propagator |T5|. It is almost equal to the usual 
lattice one except at zero and one lattice distances where the finite monopole size 
gives a modification. The first can be absorbed into a redefinition of the 'mass' term 
mo. The existence of the finite monopole size leads to two different -D(l) when the 
two currents are separated on the same line or on the parallel line. This may explain 



Note that b represents a scale at which behaviors of the monopole currents are examined. One 
should not confuse it with the physical monopole size Rm- 
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the data /2 7^ /a [0. The sohd hnes in Fig. Fig. |] and Fig. ^ are the prediction 



given by the action with the parameters written in the figures. The existence of the 
'monopole mass' mo and the running couphng constant g{b) is characterisitic of the 
action of SU(2) QCD in comparison with that of compact QED. 



7. The action if actually correct, depends only on b, not on the extendedness of the 
monopole currents. As shown above, considering the extended monopoles corresponds 
to making a block spin transformation on a dual lattice. Hence the action ( p9D is 
independent of the block spin transformation. This is similar conceptually to the 
perfect lattice action by Hasenfratz and Niedermayer ||21[ . 



8. It is interesting to see the relation between the monopole dynamics and the decon- 



finement transition in a finite volume [^. This will be discussed in a separate paper 



23 . 



In summary, the monopole condensation is shown, in the maximally abelian gauge, 
to occur always in the infinite- volume limit. The abelian charge is confined due to the 
dual Meissner effect. The U(l) charge confinement after abelian projection means color 
confinement in SU(2) QCD as proved in [Q. 

However we have not here clarified the origin of (3c due to the finite-size effect in terms 



of abelian monopoles. This will be published elsewhere [^. In addition, we have not 
yet understood what happens in other abelian projections. If abelian monopoles play an 
essential role in the continuum, the mechanism of monopole condensation may not depend 
on a gauge choice. This will be studied in future. 

We wish to acknowledge Yoshimi Matsubara for useful discussions especially on the 
monopole dynamics and the finite-size effects. Also we are thankful to Osamu Miyamura 
and Shinji Hioki for discussions on the perfect action on a lattice. This work is financially 
supported by JSPS Grant-in Aid for Scientific Research (B) (No. 06452028). 
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TABLES 

TABLE I. The quadratic terms of the monopole action adopted. All terms in which the 

relation of the two currents is equivalent are added to make each Si[k] invariant under translation 
and rotation. Here a = i> + p, b = O + p + lo, c = /i + z> + /5 + a), d = 2p, + where jl, i>, p and a) 
denote unit vectors in four different directions. 



i Si[k] 

1 kn{s)k^{s) 

2 ki^{s)kn{s + p) 

3 kfj,{s)k^{s + i)) 

4 k^{s)kfj,{s + p + !>) 

5 ki^{s)k,j,{s + a) 

6 k^{s)k^{s + 2p) 

7 kij,{s)k^{s + p + d) 

8 k^{s)ki,{s + b) 

9 k^is)k^{s + c) 

10 k^{s)k^{s + d) 

11 k^{s)k,,{s + d) 

12 k^{s)ki,{s + d + p) 
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FIGURES 




FIG. 1. An extended cube on which an 2^ extended monopole is defined as the sum of eight 
(=2^) smallest monopoles. 



15 



2z 



1- 







1 

- 


1 


1 


1 


1 1 

In7 
■III 

/ 


' 1 


e 


1 

o 










• 


o 


8 




- 

• 

o 


• 

o 


• 

o 


• 

o 


o 


□ f2 




- 
- 












^4 






♦ 

- ^ 


* 




* 


« 


9 




□ 

❖ - 


1 


1 


♦ 

1 


1 


1 1 


t 

1 


1 1 


1 



BETA 



FIG. 2. Coupling constants /j versus f3 in the strong coupling region. White symbols denote 
the SU (2) data, whereas black ones show the data of the Wilson form of compact QED with the 
substitution = (3/2. 
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FIG. 3. Coupling constants fi versus lattice size . 
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FIG. 4. Coupling constants /i versus /3 for 2^,3^, and 4^ extended monopoles on 24^ lattice. 
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FIG. 5. Coupling constants /i versus b. 
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FIG. 6. Coupling constants /2 versus b. The data with too large errors are omitted. 
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FIG. 7. Coupling constants /a versus b. The data with too large errors are omitted. 
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